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Abstract. A major challenge of today's quantum communication systems lies in 
the transmission of quantum information with high rates over long distances in the 
presence of unavoidable losses. Thereby the achievable quantum communication 
rate is fundamentally limited by the amount of energy that can be transmitted per 
use of the channel. It is hence vital to develop quantum communication protocols 
that encode quantum information as energy efficiently as possible. To this aim we 
investigate continuous- variable quantum teleportation as a method of distributing 
quantum information. We explore the possibility to encode information on 
multiple optical modes and derive upper and lower bounds on the achievable 
quantum channel capacities. This analysis enables us to benchmark single- mode 
versus multi-mode entanglement resources. Our research reveals that multiplexing 
does not only feature an enhanced energy efficiency, leading to an exponential 
increase in the achievable quantum communication rates in comparison to single- 
mode coding, but also yields an improved loss resilience. However, as reliable 
quantum information transfer is only achieved for entanglement values above a 
certain threshold, a careful optimization of the number of coding modes is needed 
to obtain the optimal quantum channel capacity. 
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1. Introduction 

Quantum communication refers to the process of transferring quantum information 
between two parties commonly called Alice and Bob. This information transfer 
forms the cornerstone of many quantum information technologies, most importantly 
quantum cryptography [I] [2] , enabling secure communication, quantum dense coding 
[3], boosting the data rates with respect to classical transmission and quantum 
networking 0]. A major challenge in all these quantum communication protocols 
is to achieve high rates over long distances in the presence of unavoidable losses. For 
this purpose, we investigate continuous- variable (CV) quantum teleportation 0(5], as 
an established method of transferring an unknown quantum state between two parties, 
using only entanglement and classical communication, which was originally introduced 
in 1993 by Bennett et al. [7] in the discrete variable regime. 

In general, all quantum communication protocols are limited by the amount of 
energy that can be transferred between the sender (Alice) and the receiver (Bob) 
per use of the channel. Consequently, the challenge in quantum communication 
resides in encoding the information as energy efficiently as possible without sacrificing 
loss resilience. For this purpose, we expand the standard single-mode CV quantum 
teleportation protocol to incorporate multiplexing. Our research shows that by 
encoding the information on multiple instead of a single mode the information transfer 
is not only more energy efficient, leading to exponentially enhanced quantum channel 
capacity in comparison to the standard single-mode protocol, but it also features an 
enhanced loss resilience. 

Furthermore, we propose a practical setup to implement the proposed 
multiplexing by encoding the information on ultrafast optical pulse modc4il. There 
exists a wide variety of sources capable to create the required entangled states 
suitable for CV quantum teleportation, ranging from optical parametric oscillators 
[H [9j [10] over four-wave-mixing in optical fibers featuring a nonlinearity [TTJ [12] 
to parametric down-conversion (PDC) in nonlinear X (2) crystals [HI HH HS1 HSl [T^ . We 
employ — without loss of generality — an ultrafast pumped PDC source which creates 
a set of Einstein-Podolsky-Rosen (EPR) states into ultrafast orthogonal frequency 
pulse modes, that can directly be applied for multiplexed quantum teleportation. 

We structured this paper into three main parts. In sections [2] and [3] we 
review the standard single-mode CV quantum teleportation protocol to introduce all 
necessary concepts and formulas. Section [4] extends the standard protocol to include 
multiplexing. In section [5] we compare the achievable quantum communication rates 
in the multiplexed regime with the standard single-mode teleportation. Section [B] 
concludes the paper and summarizes our findings. 

2. Single-mode CV quantum communication 

Before we present our multiplexed quantum communication protocol we first briefly 
review the established single-mode CV quantum teleportation scheme and the 
corresponding achievable quantum communication rates in order to introduce the 
required concepts and formulas. 

X Ultrafast optical pulses are extremely short light pulses featuring durations in the femtosecond 
regime. Using these as carriers of quantum information enables the rapid succession of states in the 
transmission further boosting the quantum communication rate. 
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2.1. Teleportation as a quantum channel 

The standard single-mode CV quantum teleportation protocol [5J [5] is illustrated in 
figure [TJ Alice intends to teleport a (unknown) quantum state p from her side to Bob. 




Figure 1. Sketch of the standard single- mode CV teleportation protocol. An 
EPR state in conjunction with classical communication is used to transmit an 
unknown quantum state from Alice to Bob. 



To this aim, Alice and Bob share a bipartite entangled state — in most cases a finitely 
squeezed EPR state — associated with the operators {a, a 1 '} on Alice's side, and 

1 6, on Bob's side obeying canonical commutation relations [d,a^] = b,& = 1. 
We denote the corresponding conjugate quadrature operators by qA = (a + a^) /v2, 
PA = (a- a f ) ji\f2 and q B = (b + 6M /v% Pb = (b- /iV2 for Alice and Bob, 
respectively. 

The CV teleportation protocol works as follows: Alice first superimposes her 
part of the shared bipartite state — we label it x — with the to be teleported state 
p. She then measures the resulting quantum system on her side and transmits the 
measurement result through classical communication to Bob. According to the data 
retrieved from Alice, Bob subsequently performs local operations on his part of the 
bipartite state x an( i obtains the teleported state ptei- 

In the scope of this paper, we are not interested in the details of the apparatus; 
hence we regard the whole protocol as a quantum channel which enables us to send a 
(unknown) quantum state p from Alice to Bob. Then, we characterize the quantum 
channel defined by the teleportation protocol in terms of its quantum communication 
capacity. A reformulation of CV quantum teleportation as a quantum channel has 
been introduced by Ban et al. [18], extending that of Bowen and Bose [19] on qubit 
teleportation. According to [18] the CV teleportation protocol with arbitrary resources 
is formally described as a generalized thermalizing channel 4>{p) = ptei, in which 
thermal-like noise decreases the teleportation qualitj{s| 

<j>(p)= dxdyf(x,y)V(x,y)pV 1 '(x,y) (1) 



§ For the qubit teleportation channel, the use of non ideal resources induces depolarization |19| . 
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Here T>(x,y) is the displacement operator 

T)(x,y)(q + ip)&(x,y) = {q - x) + i(p - y), (2) 

which shifts the input state p in its quadratures q and p according to the function 
/(x, y) given by the structure of the channel. Consequently, Bob will receive the input 
state from Alice plus some extra phase-space displacements depending on the exact 
form of the CV teleportation. The input state from Alice is distorted from its original 
form. The exact structure of the mapping function f(x, y) is dependent on the shared 
bipartite state x and is defined as 

f(x,y) = Ti{[l®'D(x,vj\ (\EPR*)(EPR*\) [i <g> ©t ( s> „)] , (3) 

where \EPR*) denotes the not- normalized EPR state 

\EPR*) = {2n)- 1 ' 2 T dq\q) A \q) B , (4) 
J — 00 

and \q) A , \q)s are the eigenstates of the quadrature operators, q A \q) A = q\q)A, 

qB\q)B = q\q)B- 

Perfect teleportation is achieved for an infinitely squeezed EPR state x = 
\EPR)(EPR\, which yields f{x 7 y) = 5(x)S(y). Hence, the input state p from Alice is 
transmitted to Bob with unit fidelity, <f>(p) = p. 

2.2. CV teleportation with Gaussian resources 

In the remainder of this paper, we restrict ourselves to a Gaussian resource x shared 
between Alice and Bob, as it is the case for the most common CV entangled state, the 
EPR state. The Gaussian state x is conveniently described in the Wigner function 
representation: 

1 

(2?r) Vdet^ 

= G(m,y)(qA,PA,QB,PB), (5) 

where £ is defined as the vector £ = (qA,PA,qB,PB), rn labels the first-order moments 
and 7 the second-order moments or covariance matrix (CM) of the state x, which 
completely characterize the Gaussian state. We have introduced the short-hand 
notation Gf m> -y) in ©, where G marks the function as Gaussian in its variables, and 
the subscripts m and 7 inside the brackets identify the first and second-order moments 
of the state. 

The first step towards evaluating the output state of the teleportation channel is to 
derive the explicit form of the noise function /(a;, y) for a given Gaussian teleportation 
resource x- Starting from the general form of f(x,y) in ([3]) the function is given by 
the convolution integral 

f(x,y)=TT J d£W E PR*{qA,PA,qB - X,p B - y)G(m,y){qA,PA,qB,PB), (6) 

where g?£ = dqAdpAolqB dps and Wepr* denotes the Wigner function of the not- 
normalized EPR state in 

To compute the convolution integral in ([6]), it is convenient to change to the 
collective quadratures ((?_,£>_, g_|_,p + ), defined as 

qA ±qB PA± PB fr7 ~. 



W x (q A ,PA,qB,PB) = /jTTT- eX P 
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In terms of the collective variables, the Wigner function of the teleportation resource 
X now reads 

W x (q-,p-,q + ,p + ) = G^- ) (q-,p-,q + ,p + ) , (8) 
where rh = (m q _ , m p _ , m q+ , m p+ ) = mR, with 
If l 2 1 2 



la being the unit matrix of size 2, and 

Iq-q- %-p- %-q+ 7q-p+ 



7 = i? J 7i? 



lp-q- 1p-p- lp-q+ 7p-p+ 
7g+g- Iq+p- lq+q+ lq+p+ 



(10) 



V 7p+9- 7p+p- 7p+<?+ 7p+p+ / 
In terms of the collective variables, the Wigner function of the not-normalized EPR 
state in dU reads 

M / BPi?*(g-,P-,'7+,P+) = 2n5{q-)6(p + ) . (11) 
We arrive at the final form of the mapping function f(x,y) for shared Gaussian 



resources 



f(x,y) = ~ [ d£5(q- +x/V2)8(p+ - y/V2)G(Thtf)(q-,p-,q + ,p + ) 



= ^G(m f ,y f )(x/V2,y/V2) 

= G (v / 2m / ,2 7/ )( a; 'y)' ( 12 ) 

where m/ = {fh q _ , m p+ ) and 

7/=(^ 9 - ) • (13) 

V ip+q- ip+p+ / 

This gives us a convenient closed formula for f{x, y) defined by the first moments m 
and CM 7 of the shared resource x between Alice and Bob. In particular, given a 
Gaussian state p on Alice's side with Wigner function 

Wp(q,p) = G {mpt ^ ) (q,p), (14) 

the teleported state </>(p) arriving at Bob's side evaluates to 



w <t>(p) P)= I dx dy f(x, y) G (mpi7p) (q-x,p-y) 

dx dy G (V2 - m/:27/) (x, v)G {mp . lp) (q-x,p-y) 

= G (m p + v / 2m / , 7p +2 7/ )( l ?'-P) • (15) 

Equation (TT5|) fully determines the CV teleportation process in the Gaussian 
framework (i.e. teleportation of Gaussian states using Gaussian resources). The 
transformation of the Gaussian input state through the teleportation channel can be 
calculated by adding the first moments and CM of the channel to the first moments 
and CM of the initial state. In the limiting case of a perfect teleportation both y/2m f 
and 27/ are zero and the initial state is retrieved. 
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2. 3. Information theoretical characterization of CV quantum teleportation 

There exist different figures of merit to quantify the accuracy of the CV teleportation. 
Among others there is the fidelity of the quantum teleportation, detailing how closely 
the state arriving at Bob's side resembles the original state from Alice. Another 
example is the classical communication capacity, giving the amount of classical 
information that can be pushed through the teleportation channel. In general, the 
choice of a figure of merit is motivated by its operational meaning. 

In the scope of this paper, we characterize the teleportation channel in terms of its 
quantum capacity [20] [21] , this means the highest rate at which quantum information 
can be reliably transmitted through the channel when Alice and Bob make use of 
error correction to convey quantum information through the noisy channel. This 
choice seems to be the most natural and appropriate, if quantum teleportation should 
be used to establish a true quantum link. 

For comparison purposes, we consider the two-way distillable entanglement as 
another figure of merit in | Appendix B| In this scenario, Alice and Bob also exchange 
classical information in a two-way fashion to extract maximally entangled states. In 
the main part of the paper, however, we will not allow two-way classical communication 
between Alice and Bob, because this approach delivers tighter bounds on the properties 
of the required resources. 

Indeed, the thermal-like noise added by the non-ideal teleportation can be 
counteracted by employing quantum error correction codes. These can increase the 
quality of the communication (e.g. in terms of the fidelity) at the cost of reducing 
the communication rate. The highest rate of reliable quantum communication, i.e. 
allowing asymptotically unit fidelity, is by definition the quantum capacity of the 
teleportation channel. The quantum capacity of Gaussian channels has been widely 
studied and characterized from an information theoretical perspective [22] [23] . In full 
generality, the quantum capacity of a quantum channel <f> is given by the following 
expression [2T] 120]: 



where <f>® indicates I parallel uses of the quantum channel. The entropic function 



is known as the coherent information. Here, S denotes the von Neumann entropy, 
S[p] = — Tr(plnp) (measured in g-nafjjj). \ip)p is a purification of p, involving an 
auxiliary quantum system denoted C, and idc is the identity quantum channel acting 
on C . In general, it is very hard to evaluate the quantum capacity of a given channel, 
because one has to optimize (|17[) over all possible input states p in the limit of infinite 
uses of the channel cf>. An analytic formula for the quantum capacity is only known 
for few specimens of CV quantum channels |23| . It is however possible to evaluate 
upper and lower bounds of the quantum channel capacity. 
In the following we put 




(16) 



i(^ e ,p) = s[^ e (p)]-s[(^ e ^id c )mM)} 



(17) 




(18) 



|| In order to obtain compact formulas for the quantum channel capacity bounds, we use natural 
logarithms, In = log e . 



Exponentially enhanced QC rate by multiplexing CV quantum teleportation 7 



This thermal-like form for the channel CM is the relevant one in several cases, as 
for the finitely squeezed EPR states with and without losses, where the parameter N 
contains the entanglement properties of the resource state. 

2.3.1. Lower Bound A lower bound on the quantum capacity can be obtained by 
restricting ourselves in (|17p to maximizing over Gaussian states pa, and by considering 
only a "single-use" of the channel, i.e. 

Q > max jo, sup 7(0, p G ) | =: Q G . (19) 

Clearly, a lower bound on the quantum capacity still provides an achievable rate of 
reliable communication^]. This lower bound can be computed efficiently for Gaussian 
channels [35] . For the teleportation channel, it is a function of the noise CM in (TT3)) . 
For a thermal- like noise with CM (fT5|) . such a quantity was computed in [32], yielding: 

Qa = max{0, —1 — IniV} . (20) 
The derivation of ([20]) is presented in |Appcndix A| 

2.3.2. Upper Bound An upper bound on the quantum capacity can be calculated 
by noting that the thermal-like noise with CM (p~8|) . for N ^ 1, can be simulated by 
the action of a linear amplifier with amplification factors I/77, followed by a linear 
attenuating channel with attenuating factor rj. In fact, the composition of these 
channels transforms the input CM 7 P to 

which coincides with the thermal-like channel by setting n = 1 — N. Due to the fact 
that the composition of channels cannot increase the quantum capacity, the capacity 
of the thermal-like channel is upper bounded by that of the attenuating channel. 
Using the results of [22] we can write 

Q < max{0, In (1-JV)- In iV}=:QA- (22) 
3. Single-mode quantum channel capacity analysis 

With formulas ([2"0]) and (|22p . we are now able to evaluate bounds on the available 
quantum channel capacities of the standard one-mode quantum teleportation protocol. 

At first, we assume that the shared bipartite entangled state is a finitely squeezed 
EPR state, 

W)pdo = ex P [ r («^ f ~ 1°) > ( 23 ) 
where the parameter r describes the generated squeezing amplitude (we assume 
without loss of generality r ^ 0), which can be transformed into the squeezing value by 
the relation: squeezing [dB] = — 101og 10 (e~ 2r ). Secondly, we study the effect of losses 
in the quantum capacity of the teleportation channel by assuming that the modes 
{a, a^}, 1 6, are attenuated by a factor rj. 

% For the case of Gaussian channel, a natural conjecture is that Gaussian states saturate the 
maximization in 1161 . However, it is in principle possible that the coherent information has a global 
maximum on non-Gaussian states. Moreover, as the coherent information might be super-additive 
for parallel channels, the regularized limit over n is in general necessary for computing the quantum 
capacity 1241 . 
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3.1. Quantum channel capacity without losses 

If we neglect losses, which can occur during the EPR state distribution to Alice and 
Bob, the parameter in the CM (fT8")l reads N = e~ 2r , where r labels the squeezing 
amplitudes of the shared EPR state. The bounds on the quantum channel capacities 
in ([20| and (|22|) evaluate to the expressions: 

Q G = max{0, 2r - 1} , (24) 

Q A = max{0, 2r + In (1 - e~ 2r )} . (25) 

The limiting factor in the CV teleportation protocol is that EPR sources are 
constrained by the maximum amount of entanglement, and hence energy, that they 
are able to emit. For the case of PDC processes, this is equivalent to the overall optical 
gain of the down-conversion process. Furthermore, the channels used to transmit the 
EPR states to Alice and Bob are constrained by the amount of energy, that they can 
carry. For example, in the case of the ubiquitous optical fibers, the most prevalent 
method for quantum state distribution, transmitted pulses exceeding a certain power 
level undergo nonlinear optical processes in the fiber and subsequently lose part of 
their entanglement. 

It is hence vital to develop quantum communication protocols that encode 
quantum information as energy efficiently as possible. For this purpose, we benchmark 
quantum communication by evaluating the quantum channel capacity as a function of 
the energy, i.e. mean photon number (n p h) inside the channel. In the case of an EPR 
state this mean photon number is given as 

(nph) =sinh 2 (r). (26) 

Figure [2] displays the calculated upper and lower bounds Qa and Qg, as defined in 
(1241) and (|25j) as a function of the mean photon number (n p h) inside of the channel. 




Figure 2. Upper Qa and lower Qq bounds for the quantum channel capacity 
(measured in q-nats) of CV quantum teleportation using a single-mode EPR 
state. The minimum squeezing required in order to reliably transmit quantum 
information resides between 3.01 dB and 4.34 dB. (K=l: A single EPR state is 
transmitted.) 

This figure shows the minimum requirements for an EPR state to enable reliable 
quantum information transfer of the teleportation channel. The upper bound Qa 
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remains zero up to mean photon numbers {n p h) — 0.125 corresponding to squeezing 
values of 3.01 dB, whereas the lower bound Qq is zero up to (n p h) ~ 0.27 equivalent 
to 4.34 dB of squeezing. Hence the minimum squeezing in EPR state allowing reliable 
quantum information transfer resides in the range between 3.01 dB and 4.34 dB. 
The situation changes if additional resources — like unbounded two-way classical 
communication — are allowed (see discussion in |Appcndix B[ ) . 

3.2. Quantum channel capacity including losses 

Analyzing quantum teleportation in the framework of realistic applications, for 
example the ubiquitous quantum state x distribution through optical fibers, the 
impact of losses has to be considered. We model these losses by the standard 
beam splitter interactions, a — > y/fja + y/1 — r]v a , b — > y/fjb + ^/l — rjiib during the 
distribution of the state to Alice and Bob, as displayed in figure [3j and evaluate the 
robustness of the state distribution as a function of the transmissivity of the channel 
r\. With these conditions N = rje~ 2r + (1 — rj), and the channel capacity formulas 




Figure 3. CV teleportation setup including standard beam-splitter like losses 
during the distribution of the EPR state to Alice and Bob. 

evaluate to: 

Q G = max{0, -1 - In [1 - ?/(l - e - 2r )}} , (27) 
Qa = max{0, In fo(l - e~ 2r )] - In [1 - 77(1 - e~ 2r )]} . (28) 

Figure [4] depicts the quantum channel capacity as a function of the transmissivity rj 
for an EPR state with a mean photon number of {n p h) = 30. 

Starting from a quantum channel capacity between 4 and 5 q-nats it quickly 
degrades for decreasing transmissivities r\ until it reaches at loss rates exceeding 
50%. 

4. Multi-mode EPR state generation and teleportation 

Having reviewed and established CV teleportation and the corresponding quantum 
communication rates in the single-mode regime we now expand the protocol to 
incorporate multiplexing. 
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Figure 4. Upper Qa and lower Qq bounds for the quantum channel capacity 
(measured in q-nats) as a function of the transmissivity r) for CV quantum 
teleportation using a single-mode EPR state including loss. The quantum channel 
capacity quickly degrades under loss until it reaches zero at loss rates exceeding 
50%. (K=l: A single EPR state is transmitted.) 

As discussed in the introduction there exist a variety of sources to create multi- 
mode EPR states. In the scope of this paper, we will focus on the properties of 
parametric down-conversion as a source of pulsed multi-mode EPR states in ultrafast 
frequency modes [551 ■ Yet our findings could also be adapted to other methods of 
squeezer generation as well. 




Figure 5. (a) Schematic PDC process: An incoming ultrafast pump pulse is 
down-converted into two squeezed and entangled signal and idler waves, (b) 
The generated PDC state incorporates a multitude of EPR states in orthogonal 
ultrafast pulse modes. 

Figure [5ja) sketches the state generation process. An incoming ultrafast pump 
pulse decays inside a medium with a x^-nonlinearity into two beams usually labelled 
signal and idler, which represent the two modes of the generated finitely squeezed 
EPR state. These states are well suited for quantum teleportation as they enable high 
repetition rates due to the ultrafast nature of the created pulses. 

However, this PDC process pumped by a pulsed laser system produces not only 
a single EPR state but, as sketched in figure [Sf a) and (b), a multitude of ultrafast 
finitely squeezed EPR states into broadband frequency pulse modes. Each output 
pulse consists of a multitude of EPR states in different orthogonal modes [27l [28], 
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formally described as 



I PDC 



lexp 



fe=i 



r k ( A\B{ - A k B k 



|0) 



(29) 



where and B k label the different ultrafast pulse modes in the signal and idler arms, 
and the parameters r k > describe the generated squeezing amplitudes. A detailed 
derivation of (|29[) is given in |25| . For common PDC sources the squeezing parameters 
r k form an exponentially decaying distribution, which can be engineered from emitting 
a single EPR state to creating a whole array of twin-beam squeezed states (see [2"9"]). 





Figure 6. Performing quantum teleportation using multi-mode PDC states in 
conjunction with multi-mode detection and displacements on Alice and Bob's side, 
effectively multiplexes the teleportation protocol. 

The standard protocol for single-mode CV teleportation [B] requires CV Bcll- 
mcasuremcnts, one-way classical communication and local phase-space displacements. 
In order to multiplex the teleportation protocol, these operations have to be performed 
on several pulse modes in parallel. There is a certain arbitrariness in that, because 
in principle different multi-mode orthogonal basis sets can be chosen for the quantum 
information encoding by the communicating parties Alice and Bob. However, in the 
following we are focusing on broadband entangled states produced via PDC, for which 
a unique natural mode basis A k B k arises from the Schmidt decomposition as given in 
(1291) . In this basis each pair of modes A k and B k forms a finitely squeezed EPR state 
and we hence can treat each teleportation independently from the others. One could 
in principle also perform the teleportation in a different basis; this however would lead 
to correlations between all individual modes, reduce the individual mode entanglement 
and consequently lower the overall quality of the teleportation. It is hence natural to 
conjecture that the basis of the Schmidt modes optimizes the teleportation capacity. 
A detailed discussion of this issue will be presented elsewhere [30] . 

These multi-mode PDC states arc hence optimally suited to multiplex CV 
quantum teleportation as a single source is sufficient for creating many EPR states 
in multiple orthogonal ultrafast frequency modes. The general multiplexed protocol 
is depicted in figure [6] From the source a multitude of EPR states is transmitted 
to Alice and Bob. Alice now encodes the state she wants to teleport in the {A k } 
modes of the source, superimposes the two beams at a beam splitter and then 
measures all optical modes separately. This can be implemented by cither splitting the 
frequency modes into different spatial modes [22 [3JH [331 131] and guiding the light to 
independent measurement setups or by performing multi-mode homodyne detection 
|35[ 136] . These measurement results are then transmitted to Bob who performs the 
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according displacements on each individual Bk mode. He then retrieves the teleported 
multi-mode state p te i- 

The experimental implementation of multi-mode teleportation represents the 
main challenge for a deployment of our multi-mode coding protocol. Alice has to 
implement homodyne measurements in multiple orthogonal modes simultaneously on 
exactly the same basis as imposed by the multi-mode EPR source. Furthermore, 
the phase reference of the local oscillator beams has to be kept stable over all optical 
modes. Any errors in the measurement basis or phase mismatch between the individual 
modes will decrease the quantum communication rate. The same reasoning also 
applies to Bob who has to perform phase-locked displacements in the exact same 
basis. Although experimentally challenging, this problem is already addressed by 
various researchers working on multi-mode homodyne detection |351 I36j and quantum 
pulse gates [3TJ GH 133 E3 ■ 

Eventually this approach of expanding the EPR-source and the detection 
apparatus to incorporate multiple modes allows us to perform multiplexed quantum 
teleportation. This in turn leads to several independent CV teleportation protocols 
being performed simultaneously. 

5. Multiplexed quantum channel capacity analysis 

In this section, we characterize the multiplexed CV teleportation channel in terms of 
its quantum capacity. 

We consider two remarkable settings. Firstly, we assume that the teleportation 
resource is given by the multi-mode EPR state in equation (|29|) . Secondly, we 
introduce a loss model in which each Schmidt mode is independently (and identically) 
attenuated by a standard beam splitter interaction with attenuation parameter 77. 

In both cases, the resulting multi-mode teleportation channel coincides with n 
parallel single-mode telcportations. Hence, proceeding as in section l2\3.1l and l2.3.2l wc 
obtain the lower bound on the multiplexed quantum channel capacity 

n 

Q G = ^max{0,-l- In N k } , (30) 

k=l 

and the upper bound on the multiplexed quantum channel capacity 

n 

Q A = ^max{0,ln(l-jV fe )-ln2V fc }, (31) 

k=l 

for suitable parameters Nk > 0. 
5.1. Multi-mode teleportation 

Neglecting losses during the EPR state distribution to Alice and Bob, the parameters 
Nk are given by Nk = e~ 2rk , where labels the individual squeezing amplitudes of 
the multi-mode squeezed state in (|29[) . The bounds on the quantum channel capacities 
in (f30|) and (|3"Tj) evaluate to the straightforward expressions: 

n 

k=l 
n 

Qa = max {°> 2rk + ln (! - e_2rfc )} • ( 33 ) 

k=l 
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The amount of energy of the multi-mode EPR state arriving at either Alice or Bob's 
side is related to the mean photon number in each arm given by: 

(n P h) =^sinh 2 (r fc ). (34) 
fe 

In analogy to the single-mode case, we analyze the teleportation channel as a function 
of the corresponding energy that now is expressed by the mean photon number (n p h) 
of all the modes involved in the teleportation protocol. 

In order to compare the standard single-mode teleportation with our multiplexed 
coding, we simulated a PDC source creating EPR states multi-mode in frequency, 
based on the source employed in |26j . The source is able to operate in various degrees 
of multi-modeness and is hence perfectly suited for comparison purposes. We designed 
it to produce three different PDC states with varying numbers of modes as presented in 
figure which shows the three normalized exponentially decaying mode distributions 
and their different weights, which we use for this analysis. These normalized mode 




k k k 



Figure 7. Three different squeezer distributions Afc normalized via Aj5 = 1 
with varying degrees of multi-modeness. Depending on the source properties 
states ranging from a single squeezer (a) up to a whole range of EPR states in 
orthogonal optical modes are generated. Here k labels the number of the generated 
finitely squeezed EPR state and Aj. its amplitude relative to the other modes. Afc 
can be converted to the actual squeezing amplitudes rj, via the overall optical 
gain B of the source: = BX^ . 

distributions can be directly converted to the corresponding EPR state distributions, 
by multiplying them with the overall optical gain B of the process = BXk (see 
[2"5] for details on the PDC source and [H]). We first simulate a purely single- mode 
source (figure Efa)), which only emits a single EPR state recreating the single- mode 
communication discussed in section [3] [5] . Figures [TJb) and (c) present states with 
rising multi-mode character, many EPR states generated in orthogonal pulse modes. 
Their effective mode numbers K = 1/Ylk ^fc E3 are K = 1,2 and 6, where it should 
be stressed that, due to the generation process, not all modes share the same squeezing, 
but the entanglement follows an exponential decay towards higher-order modes. 

Using (j32]) and (|33|) we derive the lower and upper quantum channel capacity 
bounds Qg and Qa for the different squeezer distributions presented in figure [7] The 
obtained quantities are plotted in figure [8] as function of the mean photon number or 
energy inside the channel. 

It is evident that the multiplexed teleportation relying on several less squeezed 
optical modes results in significantly higher bounds on the channel capacities with 
respect to the standard single-mode coding as soon as a certain energy threshold is 
exceeded. While the blue shaded area, which corresponds to single-mode teleportation, 
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Figure 8. Quantum channel capacity bounds (in q-nats) for multi-mode 
transmission. From bottom to top K = 1,2,6. Applying multi-mode EPR states 
for the teleportation gives a significant increase in the available quantum channel 
capacity as soon as a certain energy threshold is exceeded. This is due to the 
increased energy efficiency of multi-mode coding in conjunction with the fact that 
a minimum amount of squeezing has to be present in each optical mode to achieve 
positive quantum channel capacities (see section |3). 

with the complete energy being concentrated in a single mode, never reaches quantum 
channel capacities above 5 q-nats in the considered energy range, encoding information 
on multiple modes shows significantly higher quantum channel capacitiefl 

The underlying reason for this behaviour is the efficiency of the EPR state 
distribution. Following the discussion in |25j one finds that it is far more efficient, in 
terms of energy content, to utilize several EPR states with a low amount of squeezing 
than one EPR state with a high squeezing value. A similar effect is also observed in 
other contexts such as energy efficient entanglement creation [38| . quantum reading 
[SHUUI] and entanglement distribution [IT]. However, the fact that a certain energy is 
required to achieve a positive quantum channel capacity (see section[3]) counteracts the 
enhanced energy efficiency of multi-mode coding and consequently there exists a trade- 
off between using as many optical modes as possible for enhanced energy efficiency 
and sufficiently few optical modes to achieve positive quantum channel capacities. 

5.2. Optimal multi-mode coding 

In order to achieve the optimal quantum channel capacity one has to carefully balance 
the splitting of the energy into different modes. As discussed in section [3] the upper 
bound Qa will drop to zero as soon as the applied EPR state is below 3.01 dB. Hence, 
in order to maximize the quantum channel capacity of CV teleportation, one has to 
distribute the energy over as many EPR states as possible while the created EPR 
states still bear sufficiently high squeezing values. 

We analyzed the optimal number of modes for multiplexing that achieves maximal 
quantum channel capacities for a given amount of energy (mean photon number (n p h))- 
Our following discussion of the encoding into the optimal number of modes is split 

+ As an alternative to frequency multiplexing one could also transmit multiple weakly squeezed EPR 
states in succession instead of one strongly squeezed EPR state. Mathematically both approaches 
are equivalent. 
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into two parts: First we will elaborate on PDC sources that can be realized in a 
straightforward manner by use of existing setups, and discuss their optimal design. 
Then we turn our attention to the global optimum where the necessary squeezer 
distributions would require further engineering of the source. 



5.2.1. Common EPR-sources Given a common source of multi-mode EPR states — 
as presented in |26] — we optimize the capacities Qa and Qg over all possible effective 
mode numbers K for each mean photon number (n p h) under the restriction of a mode 
distributions r k given by the formula |29| : 



r k = By/1 - /xV , < n < 1 . (35) 
The results are depicted in figure^ Figure^a) shows the Qa and Qg bounds for the 
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Figure 9. (a) Qa and Qq channel capacities (in q-nats) for single-mode and 
optimal multi-mode coding given a common EPR-source. (b) Effective mode 
number K required for the optimal multi-mode coding. Adapted multi-mode 
codes achieve quantum channel capacities outperforming single-mode approaches. 

standard single-mode CV teleportation in comparison with the obtained optimized 
multi-mode coding. In the case of low energies both approaches yield identical rates. 
However given mean photon numbers above (n p h) ~ 0.94 (7.47 dB) and (n p h) ~ 2.40 
(10.61 dB) for Qa and Qg respectively the optimized multi-mode coding outperforms 
the single-mode approach in each bound individually. Finally, the lower bound Qg 
of the optimized multi-mode encoding surpasses the upper bound Qa of single-mode 
coding at (nph) « 5.37 (13.70 dB). 

Most importantly, however, the optimal coding bounds show a linear increase 
in channel capacity with energy, whereas the single-mode quantum capacity bounds 
exhibit a logarithmic growth for high mean photon numbers. Consequently, multi- 
mode coding enables an exponential increase of the quantum communication rate 
over single-mode coding in the presence of energy constraints. The effective mode 
number K corresponding to the optimal bounds in figure [9ta) are presented in figure 
[9jb). As the channel capacities they feature a (mostly) linear increase with energy. 

5.2.2. Optimal encoding with EPR-sources The main drawback of the currently 
available PDC sources emitting EPR states is that they feature exponentially decaying 
squeezing amplitudes r k for higher-order modes, as already depicted in figure [7l This 
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is not the optimal encoding because a certain number of squeezers will always reside 
below the bound to create positive quantum channel capacities. Hence, they do not 
contribute to the quantum communication rate while still occupying energy. 

We can negate this drawback by applying multi-mode EPR states exhibiting a 
flat distribution = r with a mode number K. Experimentally these states can 
be approximated by engineering the pump pulse and the phase-matching of the PDC 
process. This flat distribution offers the great advantage that all EPR states contribute 
to the overall channel capacity and no energy is lost in weakly squeezed modes with 
zero capacity. Indeed, it can be proven to provide the optimal distribution of the 
squeezing amplitudes, see | Appendix C| 

In the optimal case of flat mode distributions the formulas for Qq and Qa, as a 
function of the mode number K and mean photon number (n p h), evaluate to: 



Qg 



Q,K 



Qa = max < 0, K 



2 arcsinh 



2 arcsinh 



In 1 — exp 




1 



(36) 



2 arcsinh 



(nph) 
K 



(37) 



We analyze the achievable channel capacities in this optimized configuration by 
maximizing over the mode number K for given energies or mean photon numbers 
(n p h)- The results are displayed in figure [TU] Similar to the common EPR state 
distributions discussed in section 15". 2. II thev feature the advantage of showing a linear 
gain with mean photon number (n p h) instead of the logarithmic growth present in the 
single-mode coding case and hence an exponential growth in quantum communication 
rate. The achievable channel capacities surpass the quantum communication rates 
available using common EPR states as displayed in figure^ since no energy is located 
in weakly squeezed EPR states that do not contribute to the overall quantum channel 
capacity. 

Furthermore (|36p. enables us to directly assess the optimal number of modes 
K op t required to encode information for optimal capacity given a certain mean photon 
number (n p h)' 

K opt (Q G ) « 1.1133 (n ph ) , K opt (Q A ) « 2.7523 (n ph ) . (38) 

From equation (|38|l we conclude that for the optimum mode number the squeezing of 
individual modes stays fixed between 4.96 dB and 7.33 dB. Consequently using energy 
to achieve squeezing values above this threshold is actually detrimental for the overall 
quantum capacity and it is much more resourceful employing it to create EPR states 
in additional modes. 



5.3. Multi-mode analysis under loss 

We finally consider the impact of loss for multi-mode coding similar to the single 
mode case discussed in section [3J For a first analysis of the robustness under losses, 
we assume that all the modes are attenuated by the same attenuation factor 77. The 
more realistic setting of frequency depending attenuation will be considered elsewhere 
[30] . Under these conditions the channel capacity formulas evaluate to: 




Figure 10. (a) Qa and Qq channel capacities, measured in q-nats, for single- 
mode and optimal multi-mode coding given a flat mode distribution. (b) 
Effective mode number K required for the optimal multi-mode coding. Adapted 
multi-mode codes achieve quantum channel capacities outperforming single-mode 
approaches. 




PDC-state 



Figure 11. Frequency multi-mode telcportation setup including standard beam- 
splitter like losses during the distribution of the EPR states to Alice and Bob. 



Q G = ^max{0,-l-ln[l-77(l-e- 2rfc )]}, (39) 
fc=i 

n 

Qa= ^max{0,ln[?7(l-e- 2rfc )]-ln[l-?/(l-e- 2rfc )]}. (40) 
fc=i 

Using ([39]) and (|40|) we determine the loss resilience of the three exemplary states. 
We start by tuning the three test states to exhibit identical mean photon numbers 
(iT-ph) = 30 and study their behaviour under loss. Our results are visualized in figure 
[T2] where we plot the quantum channel capacity as a function of the transmissivity rj. 
Clearly an enhanced loss resilience is observed for the multi-mode coding with respect 
to the single-mode protocol which quickly degenerates under loss. The reason for this 
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Figure 12. The loss resilience of the quantum information transmission rate 
visualized for multi-mode and single-mode coding. From bottom to top K = 
1, 2,6. Multi-mode coding offers the advantage of an increased loss resilience and 
gives significant higher rates over almost the whole r\ range in comparison to the 
single-mode approach. 

advantage is well known: Strongly squeezed EPR states arc highly susceptible to loss 
whereas the encoding of information on multiple weakly squeezed states is much more 
robust against this type of noise (see, e.g., [55]). 

5.4- Optimal multi-mode coding under loss 

In a similar manner to the discussion in section [5?2] we search for the optimal number 
of modes to encode information yet including loss during the EPR state transmission. 




Figure 13. (a) Qa and Qq channel capacities (measured in q-nats) for single- 
mode and optimal multi-mode coding given a common mode distribution as a 
function of loss, (b) Effective mode number K required for the optimal multi- 
mode coding. Adapted multi-mode codes outperform single-mode approaches in 
the low-loss regime. 
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For this purpose we use an input state with mean photon number (n p h) = 30 and 
in dependence of the transmissivity r\ optimize the channel capacity over all possible 
input mode distributions. In figure [1]J] (a), we display the achievable rates using 
common squeezer distributions readily available in the lab, as already discussed in 
section [5.2.1l Figurc[13](b) depicts the effective mode numbers required to achieve the 
optimal coding. This analysis shows that in the case of losses the optimal squeezing 
values differ from the ones for lossless coding (see section 15. 3j) and the advantages 
of multiplexing are partially lost depending on the amount of loss in the system. 
In the low-loss regime the optimized multi-mode coding outperforms the standard 
single-mode approach. However in the case of high losses approaching 50% — the 
exact value depends on the initial energy or mean photon number — the single- 
mode coding surpasses our multi-mode approach. This is to be expected for the 
applied CV quantum communication protocol since it is not designed for transmission 
under extreme loss but for low-loss applications. Its optimal operational area is the 
transmission of large amounts of quantum information over short distances where it 
excels. For quantum communication over longer distances — without repeater stations 
— other quantum communication protocols are more suitable. 




Figure 14. (a) Qa and Qq channel capacities (in q-nats) for single-mode and 
optimal multi-mode coding given an optimal flat mode distribution as a function 
of loss, (b) Effective mode number K required for the optimal multi-mode coding. 
Adapted multi-mode codes outperform single-mode approaches in the low-loss 
regime. 

However, these results are still not optimal. For this purpose, we investigated the 
attainable quantum channel capacities using a flat mode distributions as discussed 
in section 15.2.21 The attainable rates are presented in figure [14] (a) and (b) . Again 
the optimized coding on flat mode distributions outperforms the single-mode coding 
in the low-loss regime and achieves higher rates than the use of common squeezer 
distributions. 

Next, we turn our attention to the quantum communication rates as a function 
of the energy for a constant loss rate. In figure ITBT a) , we plot the optimal multi- 
mode coding quantum channel capacities for a transmissivity of r\ = 0.8 for common 
squeezer distributions as a function of energy. The linear dependence of multi-mode 
quantum communication on energy for lossless coding (see section [5T2")) remains in this 
setting including losses during the state transmission. The single-mode coding also 
still features a logarithmic growth as a function of energy similar to the one observed 
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( a ) Channel capacities under loss (??=0. 8) ^ Optimal mode numbers under loss (77=0.8) 




Figure 15. (a) Qa and Qq (measured in q-nats) for single-mode and optimal 
multi-mode coding given a common mode distribution as a function of energy for 
a constant loss rate of r\ = 0.8. (b) Effective mode number K required for the 
optimal multi-mode coding. Even when considering losses the multi-mode coding 
shows an linear increase with energy, which constitutes an exponential increase 
over the logarithmic growth of the single-mode protocol. 
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Figure 16. (a) Qa and Qq (measured in q-nats for single-mode and optimal 
multi-mode coding given a flat mode distribution as a function of energy for a 
constant loss rate of n = 0.8. (b) Effective mode number K required for the 
optimal mulit-mode coding. Even when considering losses the multi-mode coding 
shows an linear increase with energy, which constitutes an exponential increase 
over the logarithmic growth of the single-mode protocol. 



for lossless state transmission. Consequently, the multi-mode protocol achieves an 
exponential increase over single-mode coding even in the presence of loss, as long as a 
certain minimum amount of energy is used in the communication. 

This effect is even more prominent when we consider optimal flat multi-mode 
EPR state distributions, as depicted in figure I16f a). It achieves higher quantum 
communication rates in comparison to the multi-mode coding on common squeezer 
distribution, while still featuring the linear growth as a function of energy as present 
in the lossless coding discussed in section 15.2.21 

However, to achieve the optimal quantum channel capacity, the squeezing values 
of the individual EPR states in the communication protocol have to be adapted to 
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the losses in the channel. Starting from the aforementioned 4.96 dB and 7.33 dB 
discussed in section 15. 2. 2 1 for lossless communication, rising amounts of EPR squeezing 
are required for optimal coding. The exact values, as a function of the transmissivity 
77, are depicted in figure [171 



( a ) Optimal EPR squeezing under loss 
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Figure 17. The optimal squeezing values, in r (a) and dB (b), for the individual 
EPR states in the multiplexed protocol adapted to the losses in the channel. 

In summary, even in the presence of loss, multi-mode coding not only gives an 
exponential increase in the observed quantum communication rate in comparison 
to single-mode coding as a function of energy, but also features an enhanced loss 
resilience. 

6. Conclusion 

In conclusion, we expanded the theory of CV quantum teleportation into the 
multi-mode domain and presented a practical approach to implement the proposed 
multiplexing protocol. We calculated upper and lower bounds on the attainable 
quantum channel capacities by encoding information on multiple optical modes. Our 
analysis reveals that multiplexing not only features an enhanced energy efficiency 
leading to an exponential increase in the achievable quantum communication rates in 
comparison to single-mode coding, but also gives an improved loss resilience. 

However, as reliable quantum information transfer is achieved only for squeezed 
modes above a certain threshold value, a careful optimization of the number of used 
coding modes is needed. 

Our findings show that EPR states with squeezing values between 3.01 and 
4.34 dB are required for having reliable quantum information transfer through the 
teleportation channel. Due to the energy constraints inside a quantum channel the 
optimum is reached when EPR states with squeezing values in the range from 4.96 
dB up to 7.33 dB arc employed. Creating squeezing above this bound is actually 
detrimental to the overall quantum communication rate. It is much more resourceful 
to invest the excess energy in creating EPR states in multiple optical modes. 
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Appendix A. Calculation of the lower bound Qc 



For computing Qg, we have to maximize the coherent information over Gaussian 
states. In this case, we can assume without loss of generality that \4>)p G 1S an EPR 
state with squeezing parameter s, shared between the subsystem A and the auxiliary 
subsystem C, described by the Wigner function, 

w m PG w(lA,PA;qc,Pc) = G {Qn A C) (q A ,PA;qc,Pc) , (A.l) 

where 
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The action of the channel transmitting the state of subsystem A from Alice to Bob, 
transforms this state to 
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(A.4) 



which is known as the Choi-Jamiolkowski (CJ) state associated with the channel. 
After tracing out the C subsystem the reduced state of subsystem B takes on the 
form 

(A.5) 



W^ p) (q B ,p B ) = G( „B)(q B ,p B ) 



(A.6) 



with 

B - I f 2A r + cosh2s 
7s ~ 2 V 2N + cosh 2s 

In order to evaluate Qc we have to determine the von Neumann entropy of the two 
states in (|A.3|) and (|A.5|) . In the case of Gaussian states this is a straightforward 
calculation, because the state is defined by its CM the von Neumann entropy is 
determined by their symplectic eigenvalues [22J E2] . Then we have 

S[4>(p)} = 9(" B - 1/2) , (A.7) 
where g(w) := (w + 1) \n(w + 1) — w\nw, and v B is the symplectic eigenvalue of the 
CM 7^ . The symplectic eigenvalue is calculated from the matrix , where il = iv-i 
is the symplectic form, with 

-1 

1 



w 2 



(Ai 
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In particular, the eigenvalues of £l-f B are ±iv B . 
Similarly, 

S[(<f> ® l){\i>)M)\ = 9(v? C - 1/2) + 9{v BC - 1/2) , (A.9) 

where v± are the symplectic eigenvalues of the CM j BC , where the ±iv BC and 
±iv BC are the eigenvalues of (fi © i})j BC . 

The resulting coherent information is an increasing function of s: 

v B = N+ i cosh 2s, (A.10) 



v BC = 1 + 2N 2 + 2iVcosh2s ± 2N\/T+ N 2 + 27Vcosh2s .(A.ll) 
In the limit of an infinitely squeezed state (s — > 00), we obtain 

v B ~N+je 2s , (A.12) 

and 



,f ^±iV. (A.13) 
Finally, after straightforward algebra, we obtain 

Qg = maxjo, lim g(v B - 1/2) - giv BC - 1/2) - g(v BC - 1/2)) 
= max{0, -1 - In N} . (A. 14) 

Appendix B. Classical communication allowed 

In the main part of the paper, we have considered a scenario in which Alice and 
Bob make use of error correction to convey quantum information through the noisy 
teleportation channel. Alternatively, if they are also allowed to exchange classical 
information in a two-way fashion, they can perform a protocol of entanglement 
purification to extract maximally entangled states up to a rate equal to the two-way 
distillable entanglement |43| . denoted D2, of the CJ state (|A.3[) . Alice and Bob can 
then use the maximally entangled states to establish a perfect teleportation channel, 
allowing reliable quantum communication up to a rate Q2 = D2 [43] - The assistance 
of two-way classical communication can in general augment the quantum capacitjQ, 
i.e., Q 2 >QM- 

We then compute the logarithmic negativity of the CJ state, denoted Qe, which 
is an upper bound for D2 |44j . To compute the logarithmic negativity, first we have to 
apply the operation of partial time reversal, denoted T, on the CJ state (|A.3j) . which 
transforms the CM (|A.4|) to 

/ 2iV + cosh2s sinh2s \ 

TV BC \ - - 2iV + cosh2s sinh2s 

( - ls >~2 sinh2s cosh 2s 

\ sinh2s cosh 2s J 

Then we compute its symplectic eigenvalues: 

d± = )-\J 2N 2 + 2N cosh 2s + cosh 4s ± (N + cosh 2s) y/ AN 2 -2 + 2 coshls" . (B.l) 

* That does not hold true for one-way classical communication 1431 . 
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The logarithmic negativity of the CJ state equals max{0, — In (2d-)}. Taking the limit 
s — y oo, after straightforward algebra, we obtain 

Q E = max{0,-lniV}. (B.2) 

Finally, generalizing this expression to the multi-mode setting, and putting N k 
we obtain 



-2r k 



Q E = 2Y j r k . (B.3) 

k=l 

Figure [Bl] shows the bounds Qq ^ Q2 ^ Qe as function of (n p h). The analysis of 
subsections 15 . Hl5~3l can be repeated for the quantity Q2 leading to similar results: The 
only qualitative difference relies on the fact that the upper bound Qe is strictly non- 
zero for all non vanishing values of the squeezing. In order to maximize this bound it 
is hence optimal to distribute the energy over as many modes as possible since there 
is no trade-off between the multi-mode structure and having zero quantum capacity 

ESI. 




Figure Bl. Upper Qe and lower Qq bounds, in q-nats, for the quantum 
channel capacity of CV quantum teleportation using a single-mode EPR state 
when classical communication between Alice and Bob is allowed. 



Appendix C. Optimal squeezing distributions 

Our aim is to optimize the squeezing distribution under energy constraint. Let us 
denote 

K 

Q:=X>(r fc ), (C.l) 
fe=i 

(K integer) the function to be optimized. We want to consider general distributions, 
including those with an infinite number of non-zero squeezers (K — > 00). To fix the 
ideas, we consider the case of lossless teleportation (the extension to the lossy case 
is straightforward). Hence, the optimization of the lower and upper bound on the 
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lossless quantum teleportation capacity is recovered by identifying the function q(r) 
with 

q G (r) = max{0, 2r - 1} , (C.2) 

or 

q A (r) = max{0, 2r + log (1 - e.- 2r )} . (C.3) 

These functions are zero if the value of r is below a certain threshold. It hence follows 
that it is sufficient to consider a finite number of squeezers corresponding to values of 
the squeezing parameters above the threshold; hence we can assume without loss of 
generality that K < oo in (|C.1|) . That also allows us to substitute the function q G , q A 
with 

q G {r):=2r-l, (C.4) 
q A {r) :=2r + \og(l-e~ 2r ). (C.5) 
In order to optimize the quantum capacity bounds under the constraint 

K 

(n ph ) = y^sinh 2 r fc , (C.6) 

k=l 

we introduce the Lagrange function 

K K 

F(ri,r 2 , ...r n ,X) = ^<?(r fc ) - A^sinh 2 r fc , (C.7) 

fe=l k=l 

with A being the Lagrange multiplier, whose value is determined by {n p h}, and q stands 
for either q G or q A . Differentiating with respect of r k we get the Lagrange equations 

dq(r k ) 



which implies 



drk 

1 dq(r k ) 
sinh(2r fe ) dr k 



Asinh(2r fc ), (C.8) 



A. (C.9) 



That means that the optimal distribution is that in which the function . ,)„ — r d ^ Tk ^ 

1 smh(2r t ) dr k 

is constant for all values of k. It hence follow that the flat distribution of the squeezing 
parameters is optimal. To check the uniqueness of the solution, we first note that 

1 dq(r k ) = dq{r(n k )) 
sinh(2r/c) dr k dn k 
where r(n k ) = arcsinh^/nfe. The Lagrange equations are then rewritten as follows: 

dg " (r(nfc)) =A. (C.ll) 
dn k 

A sufficient condition for the uniqueness of the solution is that the function q(r(n k )) 
has a given concavity as function of n k . The derivatives with respect to n k , 

dq G {r(n k )) 1 



dn k y/n k (l + n k ) ' 

dq A {r(n k )) ^arcsinh^ i 



dn k e 2arcsi„h ^ _ 1 ^ nfe ( 1+nfc ) 



(C.12) 
(C.13) 
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are both monotonically decreasing functions of n^, which proves the concavity of 
q~G( r ( n k)), and <L4(K n fc))j as functions of n^. 

In conclusion, we have proven that, for any given integer K, the flat distribution 
is the unique optimal squeezing distribution over the modes, as long as all individual 
modes feature a positive quantum channel capacity. Then, the optimal mode number 
K can be evaluated for any given (n p h), yielding the expressions presented in (|38j) . 
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